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Abstract 

In this paper, we study the dissipative dynamics of the Jaynes-Cummings 
model with phase damping in the presence of a classical homogeneous 
gravitational field. The model consists of a moving two-level atom 
simultaneously exposed to the gravitational field and a single-mode 
traveling radiation field in the presence of the phase damping. We 
present a quantum treatment of the internal and external dynamics of 
the atom based on an alternative su(2) dynamical algebraic structure. 
By making use of the super-operator technique, we obtain the solution 
of the master equation for the density operator of the quantum system, 
under the Markovian approximation. Assuming that initially the ra- 
diation field is prepared in a Glauber coherent state and the two-level 
atom is in the excited state, we investigate the influence of gravity on 
the temporal evolution of collapses and revivals of the atomic popula- 
tion inversion, atomic dipole squeezing, atomic momentum diffusion, 
photon counting statistics and quadrature squeezing of the radiation 
field in the presence of phase damping. 
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1 Introduction 



Over the last forty years many theoretical investigations have been addressed 
toward the understanding of quantum dynamics of the interacting atom-field 
system in a high-Q cavity. The interest toward this research area was mainly 
spurred by the large amount of experiments revealing the appearance of in- 
triguing features of quantum radiation-matter interaction [1]. Both theoret- 
ical and experimental activities have concentrated on trying to understand 
simple nontrivial models of quantum optics involving a single atom, regarded 
as a few effective energy levels, and one or more rear resonant modes of the 
quantized electromagnetic field. The prototype of such systems, proposed 
by Jaynes and Cummings in 1963, [2] describes a two-level atom resonantly 
interacting with a single-mode quantized field. It has proved to be a theo- 
retical laboratory of great relevance to many topics in atomic physics and 
quantum optics, as well as in the ion traps [3], cavity QED [4] and quan- 
tum information processing [5]. The Jaynes-Cummings model (JCM) also 
widely used in condensed matter physics for its relevance in spintronics [6] 
which exploits the electron-spin rather than its charge to develop a new 
generation of electronic devices [7]. When the rotating wave approximation 
(RWA) is made, the model becomes exactly solvable and its dynamical fea- 
tures can be analytically brought to light revealing remarkable properties 
[8]. The discovery of interesting aspects of the JCM as well as the devel- 
opments in cavity QED experiments involving single Rydberg atoms within 
single-mode cavities, have stimulated an intense research devoted at high- 
lighting and generalizing the original idea and physical scenario presented 
by Jaynes and Cummings. In the standard JCM, the interaction between a 
constant electric field and a stationary (motionless) two-level atom is con- 
sidered. With the development in the technologies of laser cooling and atom 
trapping the interaction between a moving atom and the field has attracted 
much attention [9-18]. 

Experimentally, atomic beams with very low velocities are generated in 
laser cooling and atomic interferometry [19]. It is obvious that for atoms 
moving with a velocity of a few millimeters or centimeters per second for 
a time period of several milliseconds or more, the influence of Earth's ac- 
celeration becomes important and cannot be neglected [20]. To get a clear 
picture of what is going to happen it may be useful to refer to the equiva- 
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lence principle. It states that the influence of a homogeneous gravitational 
field on the atom moving in a radiation field can be simulated by constant 
acceleration. This means that the following situation is physically equiva- 
lent to the atom-radiation system exposed to a gravity field: An atom is 
at rest or moving with constant velocity relative to an inertial system. The 
laboratory with the radiation field attached to it moves with constant accel- 
eration. The consequence is that the radiation field reaches the atom with 
Doppler shifted frequency. Because of the acceleration this shift changes in 
time. It acts as a time-dependent detuning. A semi-classical description of 
a two-level atom interacting with a running laser wave in a gravitational 
field has been studied [21,22]. However, the semi-classical treatment does 
not permit us to study the pure quantum effects occurring in the course of 
atom-radiation interaction. Recently, within a quantum treatment of the 
internal and external dynamics of the atom, we have presented [23] a theo- 
retical scheme based on an su(2) dynamical algebraic structure to investigate 
the influence of a classical homogeneous gravitational field on the quantum 
non-demolition measurement of atomic momentum in the dispersive JCM. 
Also, we have investigated [24] quantum statistical properties of the lossless 
Jaynes-Cummings model in the presence of a homogeneous gravitational 
field. We have found that the non-classical properties are suppressed with 
increase of the gravitational field influence. 

On the other hand, over the last two decades much attention has been 
focused on the properties of the dissipative variants of the JCM. The theo- 
retical efforts have been stimulated by experimental progress in cavity QED. 
Besides the experimental drive, there also exists a theoretical motivation to 
include relevant damping mechanism to JCM because its dynamics becomes 
more interesting. A number of authors have treated the JCM with dissipa- 
tion by the use of analytic approximations [25,26] and numerical calculations 
[27-31]. The solution in the presence of dissipation is not only of theoretical 
interest, but also important from a practical point of view since dissipation 
would be always present in any experimental realization of the model. How- 
ever, the dissipation treated in the above studies is modeled by coupling 
to an external reservoir including energy dissipation. As is well known, in 
a dissipative quantum system, the system loses energy by creating a bath 
quantum. In this kind of damping the interaction Hamiltonian between bath 
and system does not commute with the system Hamiltonian. In general this 
leads to a thermalization of the system with a certain time constant. There 
are, however other kinds of environmental coupling to the system, which do 
not involve energy exchange. In the so-called phase damping [32] the inter- 
action Hamiltonian commutes with that of system and in the dynamics only 
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the phase of system state is changed in the course of interaction. Similar to 
standard energy damping the off-diagonal elements of the density matrix in 
energy basis decay at a given rate. The phase damping can well describe 
some unaccounted decay of coherences in a single- mode micromaser [33] . It 
has also been shown that phase damping seriously reduces the fidelity of the 
received qubit in quantum computers due to the induced decoherence [34]. 
The phase damping in the JCM with one quantized field mode has been 
studied [35]. The influence of phase damping on non-classical properties of 
the multi-quanta two-mode JCM has also been studied [36]. It has been 
found that the phase damping suppresses non-classical effects of the cavity 
field in the JCM. However, all of the foregoing studies have been done only 
under the condition that the influence of the gravitational field is not taken 
into account. 

In the present contribution our main purpose is to investigate the tem- 
poral evolution of quantum statistical properties of the phase damped JCM 
in the presence of a classical homogeneous gravity field. In the JCM, when 
the atomic motion is in a propagating light wave, we consider a two-level 
atom interacting with the quantized cavity-field with phase damping in the 
presence of a homogeneous gravitational field. By solving analytically the 
master equation under the Markovian approximation, the evolving reduced 
density operator of the system is found by which the influence of the gravi- 
tational field on the dynamical behavior of the atom-radiation system in the 
presence of the phase damping is explored. In section 2, we present the mas- 
ter equation for the reduced density operator of the system under Markovian 
approximation in terms of a Hamiltonian describing the atom-radiation in- 
teraction in the presence of a gravitational field. This Hamiltonian has been 
obtained based on an su(2) dynamical algebraic structure in the interaction 
picture. In section 3 we obtain an exact solution of the JCM with the phase 
damping in the presence of a gravitational field, by which we investigate the 
dynamical evolution of the system. In section 4 we study the influence of 
the gravitational field on both the cavity-field and the atomic properties in 
the presence of phase damping. Considering the field to be initially in a 
coherent state and the two-level atom in the excited state, we explore the 
temporal evolution of the atomic inversion, atomic dipole squeezing, atomic 
momentum diffusion, photon counting statistics and quadrature squeezing 
of the radiation field. Finally, we summarize our conclusions in section 5. 
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2 Master Equation for the phase damped JCM 
in the Presence of a Gravitational Field 



The equation of motion for the density operator of the atom-radiation system 
and reservoir, p sr (t), in the Schrodinger picture is given by [37] 

= -i[H T ,p sr (t)](h = l), (1) 

where 

H T = H s + Hr + Vsr, (2) 

with the Hamiltonian of the reservoir 

H r = J2"Ah (3) 

i 

and with the Hamiltonian of the interaction between the system and reser- 
voir 

3 

V S r = H s J2Fj, (4) 



where 



,2 



A = E «a. A = E K M> ^ = h s e |s (5) 



6j and b\ are the boson annihilation and creation operators for the reservoir 
and Ki is the coupling constant. The Hamiltonian H s in (2) for the atom- 
radiation system in the presence of a classical gravity field with the atomic 
motion along the position vector x and in the RWA is given by (h = 1) 



jr - 11 

A[exp(— iq.x)a)a- + exp(iq.x)a + a], (6) 

where a and denote, respectively, the annihilation and creation operators 
of a single-mode traveling wave with frequency cj c , q is the wave vector of 
the running wave and a± denote the raising and lowering operators of the 
two- level atom with electronic levels |e), |g} and Bohr transition frequency 
uj eg . The atom-field coupling is given by the parameter A and p, x denote, 
respectively, the momentum and position operators of the atomic center of 
mass motion and g is Earth's gravitational acceleration. It has been shown 
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[23] that based on an su(2) algebraic structure, as the dynamical symmetry 
group of the model, Hamiltonian (6) takes the following form 



^ s= 2M~ M ^ + ^ + \ + X \fti(eM-iq-Z)S- (7) 
+ ex.p(iq.x)S + ), 

where the operators 

So = l(\e)(e\ - \g)(g\),S+ = a\e)(g\-^,S- = -^=|s><e|ot, (8) 

2 Vk vk 



with the following commutation relations 

[So, S±] = ±S±, [S-,S+] = -2S , (9) 

are the generators of the su(2) algebra, the operator K = a) a + |e)(e| is 
a constant of motion which represents the total number of excitations of 
the atom-radiation and A = uj eg — uj c is the detuning parameter. The 
corresponding time evolution operator for Hamiltonian (7) can be expressed 
as [23] 

u(t) = exp(iMg£t)tfu e (t)v, (10) 

where 

v = eyip(—iq.xSo),u e = exp(—iH s t). (11) 

It can be shown that the operator u e (t) satisfies an effective Schrodinger 
equation governed by an effective Hamiltonian H' s , that is 

^ = H' s u e , (12) 

where 

H' a = A(£ g)S + \Mg 2 t 2 + gSpt + A(V^_ + \fls + ) (13) 

1 q 2 - q 2 

2 u 2M u 8M 

By using the same procedure as in [24], the Hamiltonian (13) takes the 
following form in the interaction picture 

H s = u c {tfa + ^) + X - A(tf g, t)S (14) 



+ {k{t)\[ks- + k*(t)\fks + ), 
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where k(t) is an effective coupling coefficient 



^)=Aexp(|(A(Ms) + ^)), (15) 



and the operator 



-> -> 2 

h{p]t,g)=u c -{u eg + ^ + q.gt + -^), (16) 

has been introduced as the Doppler shift detuning at time t [23]. The Hamil- 
tonian (14) has the form of the Hamiltonian of the JCM, the only modifi- 
cation being the dependence of the detuning on the conjugate momentum 
and the gravitational field. Now according to ref.[23] we consider 

p sr (t) = exp(iMg.xt)i) j: (p S r)e(t)v,v = exp(-iq.xSo)- (17) 

It can be shown that the operator (p sr )e(t) satisfies an effective Schrodinger 
equation governed by an effective Hamiltonian H s = H' s + H r + V sr , that is 

.d(p sr ) e (t) 



= H s (p sr ) e (t). (18) 



By using the method given in [23] we obtain the following equation of motion 
in the interaction picture 

^Ml = -i[H s ,pi r (t)}- (19) 

The master equation for the reduced density operator of the atom-radiation 
system, /3^(t) = Tr r p^ r (t), under the Markovian approximation with ne- 
glecting the Lamb-shift term reads as [37] 

i)p ' A ' } ---HhUh^M]]. (20) 



dt 



Now, we consider 
with 



p s (t) = ulplUi, (21) 

L> 1 = T{exp(i f' H s (t')dt')}, (22) 
Jo 

where the symbol T denotes time ordering. By using (14), (20), (21) and 
(22), and following the same procedure as in ref.[37] we obtain the master 
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equation for the reduced density operator of the system under Markovian 
approximation with neglecting the lamb shift term 

^ = -i[H s ,p s (t)} -y[Hi, [Hijsit)}}, (23) 

where H s is given by (14). In Eq.(23), 7 is a parameter which depends on 
the temperature T 

7 = A^ + 2 7 rr^ ( JH|AtH|2 ), (24) 

where 

Auj' = i dujJ(oo y- v 71 , (25) 
Jo w 

and J(cj) and k(uj) are the spectral density of the reservoir and the cou- 
pling coefficient, respectively. In the derivation of the master equation we 
have assumed that the parameter T is high enough so that the Markovian 
approximation is valid. 



3 Dynamical Evolution of the Phase Damped JCM 
in the presence of Classical Gravity 

In section 2, we obtained the master equation for the reduced density op- 
erator of the atom-radiation system under Markovian approximation in the 
presence of a classical homogeneous gravitational field. In this section, We 
now start to find the exact solution for the density operator p s (t) of the 
master equation (23) with the Hamiltonian (14). For this purpose, we apply 
the approach presented in refs. [38-40]. The formal solution is given by 

p s (t) = exp(Rt) exp(St) exp(fi)p s (0), (26) 

where p s (0) is the density operator of the initial atom-field system. The 
auxiliary super-operators R, S and T are defined through their action on 
the density operator such that 

eMRt)Ps(0) = £ { -^L(H s ) k p s (0)(H s ) k , (27) 

i, J -I 

exp(St)p s (0) = exp(-iH s t)p s (0) exp(iH s t), (28) 
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exp(f t)p a (0) = exp(- 7 ( J H'J 2 t)p s (0) exp(- 7 (Fj 2 t). (29) 

We assume that initially the radiation field is in a coherent superposition 
of the Fock states, the atom is in the excited state |e), and the state vector for 
the center-of-mass degree of freedom is \ip c .m(0)) = J d 3 p(j)(p)\p). Therefore, 
the initial density operator of the atom-radiation system reads as 



Ps(0) = Pfield{0) ® Patom(O) <8> /5 c . m (0) 

where 



Pfieldify <S> Pc.m{0) 





PfieldiO) = J2J2 w niO)w m iO)\n)(m\, 



PcraiO) = I d 6 P I d 3 p'<P*ip')<Pip)\p)(p'\, 



(30) 

(31) 
(32) 



exp(-^J-)a" 



with w n i0) = ^ Py ^= ; " . The Hamiltonian (14) can be expressed as a sum 
of two terms which commute with each other, that is, 



where 



H 1 + H 2 ,[H 1 ,H 2 ]=0 
Hi = ujcicfia + -^-), 



H 2 = -A(p,g,t)S + (k(t)yjKS- + k*(t)yjKS + ). 



In the two-dimensional atomic basis we have 

Hi = uj, 



h + \ 



Ho 







A(p,g,t) 
4 



n- 5 _ 
K*(t)a 

A(p,g,t) 
4 



(33) 

(34) 
(35) 

(36) 
(37) 



Also, the square of the Hamiltonian (14) can be expressed as a sum of two 
operators, one of them is diagonal, in the form 



iH s y = Ai+A 2 , 

where 

ii = H\ + tff 

W 2 (n + i) 2 + A 2 (n + l) + (^M) 2 




(38) 



(39) 



W 2 (n-i) 2 + A 2 n+(^Mt)) 2 
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and 



A 2 = 2HiH 2 = 2lo c 



(n-i)«(i)at _(n-i)(^MM) 



(40) 



It is easily proved that L4i, A2] = 0. Taking into account the initial condition 
(30) we define the auxiliary density operator p 2 (t) as 

p 2 {t) = exp(St)exp(ft)p s (0) (41) 
= exp(-iH 2 t) exp(--fA 2 t)pi(t) exp(-^A 2 t) exp(iH 2 t), 



where the operator p\ (t) is defined by 

pi(t) = \V(t)){V(t)\®\e){e\ 

with 



(42) 



|tt(t)> = exp(- 7 t[^(n+l) 2 +A 2 (n+l)+( A( Y ,t) ) 2 ])^(0)exp(-m^t)|n). 

(43) 

(44) 
(45) 



From (36) and (39) we have, respectively 

exp(— iuj c (h + |)) 



exp(-iHit) 

exp(— j Ait) = 









exp(-iw c (n - 5)) 



(Ai)n(n,t) 

(Ai) 22 (n,t) 



where 



(Ai)u(n, t) = exp(- 7 ^ c 2 (n + \f + A 2 (n + 1) + (^M) 2 ]), (46) 



(A&nin.l) -,x,M-Vk 2 (»-^) 2 + A 2 n+(^M) 2 ])- ( 47 ) 



Also, we can write 

exp(-7^4 2 i) 

where 



ei(n, t) e 2 (h,t)a 
e%{h,t)a) ei(h,t) 



-A(p,<?,t) 1 sinh(7t v / ci(n,t)) 



ei(n,t) = cos(7t v /ci(n,t)) -a; c ( ^^ , " y )(^ + — — : 

2 x/ci(n,t) 



2 
10 



(48) 



, (49) 



e 2 (n, f) = —2co c X(n 



1 sinh( 7 yci(n- M)) 

2 y/cjjn - l,t) 



~ ,~ ,n o l, sinh( 7 t v /c 2 (n,t)) 
e 3 (n,t) = -2c c A(n--) -^—^ , 

e 4 (n,t) = cob(W5(M) - ^ c (^M)(n - V^V^M)) 



2 ^c 2 {n,t) 



(50) 
(51) 
(52) 



with 

ci(n,i) = Wg( 



*M!>)»<A + 1)(« + I)», (53 ) 



c 2 (n,t) = w c ( 



- 2 )^"/ + ^ 2 



Similarly, we can express the operator exp(— iH$t) in the two-dimensional 
atomic basis as 



exp(— iHit) 



di(h,t) d-2,(h,t)a 
ds(n,t)a^ di(n,t) 



(55) 



where 

d\(n, t) = cos(i(( 



^M) 2 +A 2 (n+l)))-( 



A(p, g, t) , sin(t((^M) 2 + A 2 (n + 1))) 

( AM M)2 + A2( . + 1)) 



S in(t((^) 2 + A 2 (fi+1))) 
^(^^) 2 + A 2 (n+l) 
S m(t((MMf + A)) 



d 2 (n, i) = — iX 



ds(n, t) = —iX 



(56) 
(57) 

(58) 



n 



* u „ A(p,g,t) , 2 2 , A(p-,ff,t ) A sin(t((^M) 2 + A 2 n)) 
ct4(n, i) = cos(i(( ) + A n)) — ( > 



Then, from (48) and (55), it follows that 
exp(— iE^t) exp(— 7yi 2 t) = 



fi(n,t) h{h,t)a 



(59) 
(60) 
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where 



fi(n,t) = ei(h,i)di(h,t) + e 2 (h,t)d 2 (h,t), 
f 2 (n,t) = e 2 (n,t)di(h,t) + e 1 (n,t)d 2 (n,t), 
f 3 (n,t) = e 3 (n,i)d 4 (n,i) + e 4 (h,t)d 3 (n,t), 
fi{h,t) = e 4 (n, t)d±(h,t) + e 3 (n, t)d 3 (n, t). 



(61) 
(62) 
(63) 
(64) 



Substituting (42) and (60) into (41), we can obtain an explicit expression 
for the operator p 2 (t) as follows 



(p 2 (t)) i ,j = \*i(t))(Vj(t)\,(i,j = 1,2), 



(65) 



with 



|*i(t)) = /i(n,t)|*(t)),|* 2 (t)) = /3(n,t)|*(t)), (66) 

where |^(t)) is given by Eq.(43). Now, we obtain the action of the operator 
exp(i?i) on the operator p 2 (t) 



k=0 



where 



which can be explicitly expressed as follows 



fjk—l ttI 



where 



«(*) 



u_ (n,i) 



( A(^t) )2+A2( . +1)) 



K*(t) 

it g k (h,t) 



( At^t) )2+A2( . +1)) 



g k + (h,t) = u k + (h,t) + 



u-(n,t), 



t(n,t)=v k + (n,t)-^^v k _(n,t), 



with 



r±(n,t) =w c (n + -)± W( j 



*M*))2 + A 2 ( a + 1))> 



(67) 



(68) 



(69) 



(70) 
(71) 



4(M = ^(r*(M)±r*(n,i)),4(M = ^(^(n,t) ± s*(n,t))> (72) 



(73) 
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J ± (n, t) = u c {n - ±) ± sj (^M^))2 + A 2 n. (74) 

Finally, by substituting (65) and (69) into (67) we obtain the exact solution 
of the master equation (23) for the phase damped JCM in the presence of a 
classical homogeneous gravity field 



Ps(t) 



Z^fc=0 fe! M ll\ t ) <Lfc=0 k\ M 12\ t ) 
Z^fc=0 fe! M 2lW 2^fc=0 fc! M 22W J 



(75) 



where 

M^(t)= (^(n,t)^ii(t)5+(n,t)+a^ fc (n,t)^ 21 (t)^(n,t) (76) 
+ g%{n,t)i> 12 {t)v'*{h,t)tf + av'^ 

M$ 2 (t)= (^(n,t)at^ n (t)a^(n,t)+5 fc (n,t)^ 21 (t)a< fc (n,t) (77) 
+ v'}(n,t)a^ 12 (t)g k _(n,t) + ^(n, t)^ 22 (t)5 fc (n, t))|^| 2 , 

M 2 *(t) = (Mf 2 (t))t = (^(n,t)atv&n(t)^(n,t) +5 fe (n,t)^ 21 (t) (78) 
xas+(n, t) + #L fc (n, 12(*)*>- («> *)a f + t)*22(*)*>- (n, t)a))\^{p)\ 2 , 

with 

= / xr4 *-(M- ( 79 ) 

v / ( A( | ^£) )2 + A2 . 

Making use of the solution given by (75), one can evaluate the mean values 
of operators of interest. In the next section we shall use it to investigate 
various dynamical properties of the phase damped JCM in the presence of 
a homogeneous gravitational field. 

4 Dynamical Properties 

In this section, we study the influence of the gravitational field on the quan- 
tum statistical properties of the atom and the quantized radiation field in 
the presence of the phase damping. 
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4a. Atomic Population Inversion 



An important quantity is the atomic population inversion which is expressed 
by the expression 

W(t) = (a 3 (t)) = Tr atom (p atom (t)a 3 (t)), (80) 

where 

Patom{t) = Trfi e l d (p s (t). (81) 

We can rewrite (80) as follows 

W{t) = Jd 3 P J2 {i\patom(t)a 3 (t)\i) = J d 3 pf;((n|®((e|p s (t)|e)-( 5 |p s (t)| 5 ))®|n)). 

i=e,g n=0 

(82) 

Therefore, by using (75) and (82) we obtain 



W(t) = Jd^(±± i -^((n\M* l ( 



L (t)|n) - (n\M^t)\n))), (83) 
where from (76) and (77) we have 

(n\M^(t)\n) =(g k + (n,t)) 2 \Mn,t)\ 2 + (h k (n+l,t)) 2 (84) 
x |V> 2 (n + 1, t)\ 2 + 2Re[g k + (n, t)h k {n + 1, i)$f(n, t)ifo(n + 1, *)], 



(n\M k 2 (t)\n) = (h k (n,t)) 2 \Mn - l,t)\ 2 + (<£(n,t)) 2 |^(n,t)| 2 (85) 
+2Re[g k (n, t)/i fc (n, t)^(n - 1, t)ih(n, t)], 

with 

h k (n,t) = ^v' k (n,t), (86) 

and 

^ i (n,t) = H* i (t)),(i = l,2), (87) 

where we have defined |vl/j(i)) in (66). 

In Fig.l we have plotted the atomic population inversion as a function 
of the scaled time At for three different values of the parameter q.g. In this 
figure and all the subsequent figures we set q = 10 7 m _1 , M = W~ 26 Kg, 
g = 9.8p, io rec = g = .5x 10 6 ^, A = 1 x 10 6 ^, A = 8.5 x 10 7 ^, 
a = 2, A = 1.8 x 10 6 ^, <f>(p) = - 7 i=exp(^) with a = 1 [21-23] and 
7 = 7 x 10 -5 — . It should be noted that the relevant time scale introduced 
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by the gravitational influence is r a = — 3= [23]. For an optical laser with 

q = 10 7 m _1 , r a is about 10 _4 s. In Fig.la we consider small gravitational 
influence in the presence of the phase damping. This means very small 
q.g, i.e., the momentum transfer from the laser beam to the atom is only 
slightly altered by the gravitational acceleration because the latter is very 
small or nearly perpendicular to the laser beam. In Figs, lb and lc we con- 
sider the gravitational influence in the presence of the phase damping for 
q.g = 0.5 x 10 7 and q.g = 1.5 x 10 7 , respectively. By comparing Figs. la, lb 
and lc we can see the influence of gravity on the time evolution of the atomic 
population inversion when there is the phase damping. As it is seen from 
Fig.la for the atomic population inversion the Rabi-like oscillations can be 
identified. With the increasing value of the parameter q.g (see Figs. lb and 
lc) the Rabi oscillations of the atomic population inversion disappear. 

4b. Atomic Dipole Squeezing 

To analyze the quantum fluctuations of the atomic dipole variables and ex- 
amine their squeezing we consider the two slowly varying Hermitian quadra- 
ture operators 

<7i = ^((T + exp(-iuj eg t) + <t_ exp(iuj eg t)), (88) 

and ^ 

<?2 = exp(-iw e9 t) - <t_ exp(iuj eg t)). (89) 

In fact <7i and &2 correspond to the dispersive and absorptive components 
of the amplitude of the atomic polarization [41], respectively. They obey 
the commutation relation [<5"i,<5"2] = §^3- Correspondingly, the Heisenberg 
uncertainty relation is 

(Aa 1 ) 2 (Aa 2 ) 2 >^\(a 3 )\ 2 , (90) 

where (A<jj) 2 = (of) — (<7j) 2 is the variance in the component = 1, 2) of 
the atomic dipole. 

The fluctuations in the component = l,2)are said to be squeezed 
(i.e., dipole squeezing)if the variance in &i satisfies the condition 

(Aa i ) 2 <±\{a 3 )\,(i = lor2). (91) 
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Since of = \ this condition may be written as 

Fi = l- 4(ai) 2 - \(a 3 )\ < 0, (t = lor2). (92) 
The expectation values of the atomic operators <r + and <r_ are given by 

(a±(t)) = Tr atom {p atom {t)a±{t)). (93) 
We can rewrite (93) as follows 

/oo 
d 3 p]T((n| ® «e|p a (t)a ± (t)|e) - (g\p s (t)a±(t)\g)) |n)). (94) 

n=0 

Therefore, by using (75) and (94) we obtain 

<*_(*)> = / d'VE E ^TT-^l^faWIn) = <*+(*)>*, (95) 

J fc=0n=0 ^ 

where from (78) we have 

(n\M k 2 (t)\n) = g k + (n, i)Y>i(n, t)(g k _ (n, t)^{n, t) + h k (n, t)^\ (n - 1, t)) (96) 
+h k (n + 1, t)^(n + 1, t)(h k (n, t)^{ (n - 1, t) + <? fe (n, t)rl>\ (n, *)). 

The time evolution of corresponding to the squeezing of <7i has been 
shown in Fig. 2 for three values of the parameter q.g in the presence of the 
phase damping. As it is seen, with the increasing value of the parameter q.g 
the dipole squeezing is completely removed. 



4c. Atomic momentum diffusion 

The next quantity we examine is the atomic momentum diffusion. As a 
consequence of the atomic momentum diffusion, the atom experiences light- 
induced forces (radiation force) during its interaction with the radiation 
field. The atomic momentum diffusion is given by 

Ap(t) = {{vit?) - <P(i)> 2 )i (97) 
Now we calculate the expectation values of the p and p 

$(t)} = Tr atom {p atom {t)p{t)), (98) 

$(t) 2 ) = Tr atom (p atom (t)p{t) 2 )- (99) 
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By using p\p) = p\p) and (75) we obtain 



(p(t)) = M 3 PEE ^pdnlM^n) - (n\M* 2 (t)\n)), (100) 

J fc=0n=0 

(p(tf) = /^EE ^/(HM^n) - (n|M 2 fe 2 (t)|n)), (101) 

fc=0n=0 

where (n|Mf L (t)|n) and {n\M2 2 (f)\ n ) are given by (84) and (85), respec- 
tively 

In Figs. 3a-3b we have plotted Ap(t) for q.g = 0, q.g = 0.5 x 10 and 
q.g = 1.5 x 10 7 , respectively. In figure 3a the Rabi-like oscillations can be 
identified, but in Figs. 3b and 3c, when the influence of the gravitational field 
increases, the Rabi oscillations disappear and the atomic momentum diffu- 
sion becomes positive and the atom does not recoil, because before absorp- 
tion of photon, the atom is deflected by the gravitational field . Moreover, 
the atom can experience larger light-induced forces during its interaction 
with the radiation field, when the gravitational field increases. 

4d. Photon Counting Statistics 

We now investigate the influence of gravity on the sub-Poissonian statistics 
of the radiation field. For this purpose, we calculate the Mandel parameter 
defined by [42] 

For Q < (Q > 0), the statistics is sub-Poissonian (super- Poissonian); 
Q = stands for Poissonian statistics. Since (n(t)) = J2'^=o n ^'( n ^) an d 
(n(t) 2 } = E^° =0 n 2 P(n, t) we have 

00 00 00 

Q(t) = ({[£n 2 P(n,t)] - [^PK^ME^Mr 1 ) " 1, (103) 
n=0 n=0 n=0 

where the probability of finding n photons in the radiation field is found to 
be 

P(n,t) = (n\pf ie i d (t)\n) = (n\Tr atorn p s (t)\n), (104) 
and by using (75) we have 

P(n,t) = d 3 ? EE { -V-((n\M^(t)\n) + <n|M* (t)|n». (105) 

J k=0n=0 K ' 
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Therefore, by using (103) and (105) we obtain 

OO - OO OO .xfc 

= ({[£ ^( M 3 fEE nF" ( < n|M " ( * )|n) + (n|M 2 fc 2 (t)|n))X!06) 

n=0 J fc=0n=0 
oo . oo oo (n f \k 

"IE »( / ^E E ^((nK^ln) + (n|M 2 fc 2 (t)|n)))] 2 } 

n=0 J fc=0n=0 
oo . oo oo (n f \k 

x[£n(/ d 3 p£ E HA ((n|M " l( * )|n) + ^1^22 Win)))]" 1 ) - 1. 

n=0 J fc=0n=0 

The numerical results for three values of the parameter q.g are shown 
in Fig. 4. As it is seen, the cavity-field exhibits alternately sub-Poissonian 
and super-Poissonian statistics when the influence of the gravitational field 
is negligible. With increasing q.g the sub-Poissonian characteristic is sup- 
pressed and the cavity-field exhibits super-Poissonian statistics. After some 
time, the Mandel parameter Q is stabilized at an asymptotic zero value; the 
larger the parameter q.g is more rapidly Q(t) reaches the asymptotic value 
zero. 

4e. Quadrature Squeezing of the Cavity-Field 

Finally, we investigate the influence of gravity on the quadrature squeez- 
ing of the radiation field. For this purpose, we introduce two slowly varying 
quadrature operators 

X\(t) = -(aexp(iu;t) + at exp(— iu>t)), (107) 

and ^ 

A 2 (i) = — (aexp(iu>t) — at exp(— iut)), (108) 

where a and at obey the commutation relation [a, at] = 1. The operators 
Xi(t) and X 2 (t) satisfy the commutation relation 

[X 1 (t),X 2 (t)]= 1 -, (109) 

which implies the Heisenberg uncertainty relation 

((Al 1 (t)) 2 )((Al 2 (t)) 2 ) > 1. (110) 

A state of the radiation field is said to be squeezed whenever 

{(AX i ) 2 )<±(i = lor2), (111) 
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where 

((AX i ) 2 ) = (X 2 )-(X i ) 2 ,(i = l,2). (112) 

The degree of squeezing can be measured by the squeezing parameter Si,(i = 
1,2) defined by 

S i (t)=4((AX i (t)) 2 )-l, (113) 

which can be expressed in terms of the annihilation and creation operators, 
a and at as follows 

5i(t) = ((a\t)) - (a(t)) 2 )exp(2^t) + ((a^ 2 (t)} - (at(t)) 2 ) (114) 

x exp(-2iujt) + 2((a t (£)a(£)) - (a t (t))(a(t))), 

and 

S 2 (t) = -«a 2 (t)) - (a(t)) 2 )exp(2^t) - {(at 2 (p,t)) - <at(i)> 2 )(115) 

x exp(-2iwt) + 2((a t (t)a(t)) - (a^t))^*))). 

Then, the condition for squeezing in the quadrature components can be 
simply written as S{(t) < 0. Now we obtain (a(t)}, (at(i)) and (at(t)a(t)) 

(a(t)) = Tr field (p field a(t)), (a 1 (t)} = Tr field (p field at (t)), (116) 

{at(t)a(t)) = Tr field {p field a\t)a(t)), (117) 
so that by using (75) we have 

(a(t)) =/^EE { -^ L ((n\M^ 1 (t)a\n) + (n\M* 2 (t)a\n)) = (a^t))*, 

fc=0n=0 

(118) 

r oo oo (n f\k 

(a(t) 2 ) = / E ^((nlMf^tja^nJ+^IAf* (t)a 2 |n>) = <at 2 (t)>*, 

J fc=0n=0 K - 

(119) 

(at(t)a(t)) =/d 3 pEE ^-"((nl^iWI") + MM 2 fe 2 (t)|n)) (120) 

fc=0n=0 

where 

(nlM^^aln) =V5Js* + (n,t)s*( n -l 1 i)^( n ,t)^(n-l 1 i) (121) 
+^/n(n + l)v+ (n + 1, t)s+(n - 1, + 1, t) 

x^*(n- l,t) + nv+(n,t)g'l(n,t)'il>i(n,t)'il>2(n,t) 
+nVn + lv'}(n + 1, t)v'*(n, t)ip 2 (n + 1, t)V4 ("-, *)> 
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= ^n(n- l)g k + (n, t)g k + (n - 2, t)^{n, t) (122) 
il>\ (n - 2, t) + ^/n(n-l)(n+l)*/*(n + 1, t)c?+(n - 2, t) 
x^ 2 (n + 1, *)^* (n -2,t) + (n- l)y/Ev'j!(n - 1, t)s+(", i) 
xV>i(n,t)V>2(n - + (n - l)\Jn(n + l)t/_ fc (n + 
x^(n - l,*)^2(n + l,t)^2(n - 1,*), 

= v^5- (n, t)9- (n - 1, i)V>2 (n, t)V>2 (« - 1, *) (123) 
+yjn{n- l)v' k {n- l,t)g k + {n,t)^ 2 {n,t) 
X7pl(n-2,t) + nv'^(n,t)g k _(n- l,t)tpi(n - ^^(rc- 1,*) 
+n\/n - lt/_ fc (n - 1, t)v' k (n, i)ipi(n - l,t)V>|(n - 2,i), 

= ^/n(n - 1)/ (n, t)/ (n - 2, t)ih (n, t) ^ (n - 2, i) (124) 
+^/n(n-l)(n-2)v!f(n - 2, (n, t)V> 2 (n, t) 
x^*(n - 3,t) +nVn 3 Tv!? (n - 2,t) 
x^i(n - 1, *)$j(n - 2, t) + (n - 1) ^(n - l)(n - 2)u'*(n, t) 
xu!?(n-2,*)^i(n- 1, t)^*(n - 3, t), 

In Fig.5 we have plotted the squeezing parameter versus the scaled 
time Xt for three values of the parameter q.g. As it is seen, the quadrature 
component X\ exhibits squeezing in the course of time evolution when the 
influence of the gravitational field is negligible. With increase of the param- 
eter q.g, the parameter Si shows damped oscillatory behaviour and there is 
no quadrature squeezing. 

5 Summary and conclusions 

In this paper we studied the dissipative dynamics of the JCM with phase 
damping in the presence of a classical homogeneous gravity field. The model 
consists of a moving two-level atom simultaneously exposed to the gravita- 
tional field and a single-mode traveling radiation field in the presence of 
the phase damping. We presented a quantum treatment of the internal 



(n\M^(t)a 2 \n) 



(n\M k 2 (t)a\n) 



(n\M% 2 {t)a 2 \n) 



20 



and external dynamics of the atom based on an alternative su(2) dynami- 
cal algebraic structure. By making use of the super-operator technique, we 
obtained an exact solution of the master equation for the density operator 
of the quantized atom-radiation system, under the Markovian approxima- 
tion. Assuming that initially the radiation field is prepared in a coherent 
state and the two-level atom is in the excited state, we investigated the in- 
fluence of gravity on the temporal evolution of collapses and revivals of the 
atomic population inversion, atomic dipole squeezing, atomic momentum 
diffusion, photon counting statistics and quadrature squeezing of the radia- 
tion field in the presence of the phase damping. The results are summarized 
as follows: 1) the Rabi-like oscillations in the atomic population inversion 
disappear, 2) the dipole squeezing decays with increase of the parameter 
q.g, 3) in the presence of the gravitational field, the atom can experiences 
larger light-induced forces during its interaction with the radiation field , 
4) with increase of q.g, the sub-Poissonian behaviour of the cavity-field is 
suppressed and it exhibits super-Poissonian statistics and after some time, 
the Mandel parameter Q(t) is stabilized at an asymptotic zero value; the 
larger the parameter q.g is more rapidly Q(t) reaches the asymptotic value 
zero, and 5) the quadrature squeezing of the cavity-field disappears. 
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FIGURE CAPTIONS: 

FIG. 1 Time evolution of the atomic population inversion versus the 
scaled time At. Here we have set q = 10 7 m _1 , 
M = W- 26 kg,g = 9.8f,u rec = .5 x 10 6 ^, 

A = 1 x 10 6 ^, A = 8.5 x 10 7 ^, if = 0, a = 2, A = 1.8 x 10 6 ^, 
7 = 7 x 1(T 5 ^, 



a) For q.g = 0. 

b) For q.g = 0.5 x 10 7 . 

c) For q.g= 1.5 x 10 7 . 

FIG. 2 Time evolution of the atomic dipole squeezing versus the scaled 
time At with the same corresponding data used in Fig.l; 



a) For q.g = 0. 

b) For q.g = 0.5 x 10 7 . 

c) For q.g = 1.5 x 10 7 . 



FIG. 3 Time evolution of the atomic momentum diffusion versus the 
scaled time At with the same corresponding data used in Fig.l; 



a) For q.g = 0. 

b) For q.g = 0.5 x 10 7 . 

c) For q.g = 1.5 x 10 7 . 



FIG. 4 Time evolution of the Mandel parameter Q(t) versus the scaled 
time At with the same corresponding data used in Fig.l; 



a) For q.g = 0. 

b) For q.g = 0.5 x 10 7 . 

c) For q.g = 1.5 x 10 7 . 



FIG. 5 Time evolution of the squeezing parameter S±(t) versus the 
scaled time At with the same corresponding data used in Fig.l; 



a) For q.g = 0. 

b) For q.g = 0.5 x 10 7 . 

c) For q.g= 1.5 x 10 7 . 
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